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I The baryon fraction in galaxy clusters is one of the most direct way to constrain fio- The 

. baryonic fraction is estimated to be in the range /b ~ (0.15 — 0.20) h'^o^^ which is several 

times higher than expected from the observed light element abundances in f2o = 1 universe, 
leading to the conclusion that we live in low density universe. In this work I will first present 

■ the various steps which lead to these results and then foccus on the distribution of the baryon 
' fraction inside clusters. I will show that while the baryon fraction in clusters follows a scaling 
, law, its theoretical expected shape does not agree with what is inferred from observations. I 

1/^ ' will show that when various factors entering into the determination of clusters gas fraction 

' are taken into account the observed baryon fraction profile matches better the predictions, 

provided that the average baryon fraction at the virial radius is lower than previously found 
/b ~ 10%. 

^ . 1 Introduction 

H , 

I The baryon fraction in clusters, /b, is assumed to be equal to the cosmic value. Combined with 

' the universal baryon density Jib as predicted from light element abundances through the theory 

of big-bang nucleosynthesis it yields to the upper limit of the cosmic mass density: 



f^o = r(r5oo)T^^ (1) 



The numerical factor T (~ 0.92) is a correction factor for possible differences between the 
baryon content in clusters and the universal one arising during cluster formation. There are 
still uncertainties on the current BEN results on the primordial D/H abundance with 0.024 < 
^50 — 0.08. As the formulae shows, this test to JIq is in principle the most direct one as 
it depends on very few reasonable assumptions: the cluster formation by collapse from a well 
mixed medium, and no segregation between the gas and the dark matter as shown by numerical 
simulations. 

Rosat observations of Coma cluster have led to the following numbers: Baryonic content : Mf, ~ 
3.4 H^q^'^W^'^Mq; Mass in stars : ~ 2.2 H^^W^^Mq; Total mass : Mtot ~ 2.2 H^^W^^Mq. 
These numbers yields to a gas fraction of /gas ~ 0.15 h^Q at Abell radius. When combined to 
current estimates of f^b ^io ^^is leads to low value of JIq (White et al. 1993). Several analyses 
of large samples of clusters have confirmed such a high baryon fraction with values ranging 
from 15 — 20%. Hower, there is large dispersion between the published values and it has even 
been argued that the baryon fraction increases with cluster temperature (mass) and with radius 
inside a cluster. This fact has often been explained as due to non-gravitational processes such 
as heating effects occuring during galaxy cluster formation. In this talk I will address the issue 



of the baryon fraction in clusters by examining critically its distribution within the cluster and 
I will show how some systematic effects can bias the baryon fraction determination. 



2 The standard method of gas fraction determination 

The main assumption of this method is that the galaxy baryon fraction reflects the universal 
one. The main component of the baryonic content is mostly in the form of hot X-ray gas. The 
standard way to derive the gas fraction proceeds from the following steps: 

• Assume a P — model to fit the observed surface brightness So{l + {6/0c)'^)~^^~^^/'^ where 9 is the 
projected angular distance to the center. The gas density profile (thus the gas mass profile) is 
then derived using the two best-fit parameters /? and 6c- The X-ray gas emission is rarely traced 
out to very large radii, so that gas fractions are measured only up to an X-ray limiting radius 
for which the signal-to-noise is good enough. Most of the time, it was necessary to extrapolate 
to the virial radius r2oo or to some other outer radius like rsoo- It is then essential to realize 
that gas mass is estimated at a radius at which the actual emission is poorly constrained, with a 
value of the order of the X-ray background or less. Moreover, the emissivity being dominated by 
the central region, a parametric fit will be rather insensitive to the outer part of the gas profile. 

• Use the hydrostatic equation to derive the total mass inside the same radius: Mtot{r) = 

G^rrip /^^x^ ^1 + ^ ' "^^^ total mass thus depends linearly on both (3 and Tx- Hence, 

if the slope of the gas density is poorly determined, it will have a drastic influence on the derived 
mass. 



3 The baryon fraction profile: Observations versus simulations 

In the case where only gravity is acting during cluster formation, the baryon (gas) fraction profile 
is expected to follow a scaling law, depending only on the contrast density S {S = p{< r)/pc 
where pc = /8ttG) . Hydro-dynamical simulations of an X-ray cluster from different groups 
(the Santa Barbara group) have shown that the gas fraction normalized to the global value, is a 
function of radius it first increases in the inner part and then tends to flatten in the outer part to 
reach the (cosmic) value see Fig. |l] . The distribution of the gas fraction within clusters has been 
widely examined by Roussel et al. (2000) (hereafter RSBOO) for a large sample of groups and 
galaxy clusters. The resulting shape is shown in (Fig. |l|). To ensure consistency I will use the 
scaling law derived gas fractions by RSBOO, in which masses are derived from the NFW universal 
dark matter profile normalized to the Ms — Tx relation from numerical simulations with two 
calibrations: Tx = 4.75(M2oo/lO^^M0)2/3keV (EMN96) and Tx = 3.81{Mi7s/W^^ Mq)'^/^ keV 
(Bryan and Norman 1998, BN98), Ms being the mass enclosed in a region with a density equal 
to (1 + 6) times the critical density. These two values of the normalization can be considered 
to be the extreme values among existing numerical simulations. BN98 normalization leading 
to virial masses nearly 40% higher. For comparison gas fractions at rsoo and r2oo from the 
literature are also plotted in (Fig. ||). The comparison of both profiles is very surprising as it 
shows that: the gas fraction profile derived from observations is in strong disagreement with the 
numerical simulations results (for a global value of 16%), its shape continuously increases and 
does not exhibit the flattening in the outer parts as seen in numerical simulations. Clearly, this 
discrepancy calls for caution when one is using the gas fraction to set upper limit on the mean 
density of the universe. 

One possibility is that processes acting during the cluster formation are not well understood and 
probably non-gravitational processes could have played an important role. But this is probably 
not the case, because in that case large departure from scaling laws is expected which is not 
observed in the data (RSBOO). 



Figure 1: The distribution of the gas fraction /gas ver- 
sus 6. The observations are from RSBOO filled circles, Et- 
tori & Fabian 1999 empty triangle and Arnaud & Evrard 
1999 crosses. Statistical uncertainties on these quantities 
are smaller than the symbols size. The squares connected 
by a dashed line correspond to the theoretical gas fraction 
for a global value of 16% (from the numerical simulations 
of the Santa Barbara group). Small squares connected by 
a line correspond to the gas fraction calculated from nu- 
merical simulations including winds by Metzler & Evrard 
(1997): for a 6 keV cluster {long dashes), a 3 keV clus- 
ter {small dashes) and without including winds {continuous 
line) with a global value of 20%. 

4 Uncertainties in gas fraction estimation 

4-1 The extrapolation effect 

In a recent study of a large sample of X-ray ROSAT images in which the data trace the emission 
up to very large radii (approximately up to the virial radius), it has been shown that the P- 
model does not provide an accurate description of the surface brightness over the whole range 
of radii. They have found that the outer slopes (0.3r2oo to r2oo) are actually steeper than the 
slope found when the /3-model is appHed to the whole cluster which actually leads to lower gas 
masses (Vikhhnin et al. 1999, VFJ99). 

4-2 The clumping effect 

In most of the studies on gas and total mass estimations, the ICM gas is assumed to be uniform, 
which is not true, as actual clusters do exhibit a certain level of density fluctuations at small and 
large scale probably due to accretion and mergers events. This clumping is assumed to be more 
pronouced in the outer parts of the clusters where the relaxation is not completely achieved. In 
the presence of clumping the gas masses are overestimated by a factor = l.lGat 6 = 500 
(C =< pI > / < pg >2) (Mathiesen et al. 1999, MEM99). 

5 Constraining Oq using the corrected gas fraction 

I will now estimate the new gas fractions taking into account the above effects. The gas masses 
are corrected from the extrapolation problem by using the VFJ99 sample (we use their results to 
compute the gas mass in the virial region) and for the clumping effect using the MEM99 factor. 
Virial masses are computed from numerical simulations as described in section 3. The result is 
plotted in (Fig. ^). From this figure, we can see that the gas fraction distribution, in particular 
the asymptotic behaviour, is now consistent with numerical simulations predictions, (although 
a slight difference persists in the inner parts at 5 < 10^ which can well be due to non-negligible 
energy injection in these regions). By matching the data point at 6 = 500, we find /gas = 
0.0875 ±0.0075 (0.108 ±0.0092) at 68% confidence level with BN98 (EMN96) normalization, the 
uncertainty being due to numerical simulations (Fig. ^). Adding the stellar mean contribution 
of 1%, our final value of the baryon fraction at the virial radius is: /b ~ 0.10 ± 0.01 (with 
BN98 calibration) and /b ~ 0.12 ±0.01 (with EMN96 calibration). Plugging this new value 
of the baryon fraction in equation (1) and using the low D/H value as recently reported by 
O'Meara et al. (2000) corresponding to Ob ~ 0.08 (this value is consistent with recent results 




on CMB fluctuations measurements from Boomerang mission 2001), we find the following value 
for r^o- ^0 = 0.8 lb 0.1 using the baryon fraction estimated with BN98 normalization. This value 
is consistent with the result from the evolution of X-ray clusters abundance (Blanchard et al. 
2000). If we use the baryon fraction obtained with EMN96 calibration, we find ~ 0.65. 
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Figure 2: The distribution of the gas fraction fg versus 
S: circles correspond to the observed fg from RSBOO 
but restricted to the region where the gas is detected 
(no extrapolation), using EMN96 {empty circles) and 
BN98 (filled circles) calibration. Empty (filled) star 
corresponds to the estimated fg using VFJ99 sample 
and corrected for the clumping with EMN96 (BN98) 
calibration. Filled squares linked by a dashed line 
correspond to the gas fraction predicted by numerical 
simulations for a global fraction of 8.75%. The upper 
dashed line corresponds to the case where the global 
fraction is 10.8% 



6 conclusion 

The baryon fraction in X-ray clusters is one of the most direct way to constrain 0,q. Previous 
estimations of the baryon fraction yield to a high value /b ~ 15% /ijq , favouring low density 
universe 0,q ~ 0.3. In my talk I have shown that the radial profil of the gas fraction as derived 
from observations differs strongly from what is found in numerical simulations. I argue this 
is not due to non-gravitational processes which may occur during the cluster formation but is 
rather due to various systematics such as the extrapolation of a /3-model and the dumpiness 
of the gas that have not been taken into account in previous estimations. Correcting the gas 
masses from these biases and using total masses calibrated from numerical simulations yield to 
a better agreement between observations and numerical simulations for a global value of the 
order of ~ 10(12)% h^Q depending on which normalization we use to compute the total mass. 

— 1/2 

This revised value is lower than previous ones and leads to JIq ~ 0.8(0.65) /i^q 
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